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Abstract. 

In this review, I discuss briefly stringent tests of Lorentz- violating quantum space-time foam 
models inspired from String/Brane theories, provided by studies of high energy Photons from 
intense celestial sources, such as Active Galactic Nuclei or Gamma Ray Bursts. The theoretical 
models predict modifications to the radiation dispersion relations, which are quadratically 
suppressed by the string mass scale, and time delays in the arrival times of photons (assumed 
to be emitted more or less simultaneously from the source), which are proportional to the 
photon energy, so that the more energetic photons arrive later. Although the astrophysics 
at the source of these energetic photons is still not understood, and such non simultaneous 
arrivals, that have been observed recently, might well be due to non simultaneous emission 
as a result of conventional physics effects, nevertheless, rather surprisingly, the observed time 
delays can also fit excellently the stringy space-time foam scenarios, provided the space-time 
defect foam is inhomogeneous. The key features of the model, that allow it to evade a plethora of 
astrophysical constraints on Lorentz violation, in sharp contrast to other field-theoretic Lorentz- 
violating models of quantum gravity, are: (i) transparency of the foam to electrons and in general 
charged matter, (ii) absence of birefringence effects and (iii) a breakdown of the local effective 
lagrangian formalism. 



1. Introduction: Mysterious Results in High Energy Gamma- Ray Astronomy 

On July 9th 2005, the MAGIC (Afajor Atmospheric Gamma-ray /maging Cherenkov) 
Telescope, located in the Canary Islands, observed pQ very high energy gamma rays from the 
active galactic nucleus Markarian 501 (Mkn 501), which lies at red-shift z = 0.034 {i.e. about 
half a million light years away) from Earth, with energies up to the order of 10 TeV (1 TeV = 
10 3 GeV = 10 12 eV), which were delayed up to four minutes as compared with their lower-energy 
counterparts (in the 0.6 TeV or lower range). It was the first observation of such a distinct delay 

Three years after the MAGIC observations, in September 2008, the FERMI (formerly known 
as GLAST) Satellite Telescope pj, also observed time delays of the higher-energy photons, from 
the distant Gamma Ray Burst (GRB) 080916c [3], at red-shifts z = 4.35, and later on from 
GRB 090510 @], at red-shift z = 0.9 and from GRB 09092B, at redhifts z = 1.822 [5]. 

The delay effects may be due to the conventional astrophysics of the active galactic nucleus 
or Gamma- Ray Burst {source effect), which, however, is not well understood at present [HE]. 
In fact, currently there seem to be no consensus among the astrophysicists on the appropriate 
mechanism for the production of such high-energy photons at the source. 



These uncertainties prompted more ambitious, although admittedly far-fetched, explanations 
for the MAGIC effect [$\ and more general for the other observed delayed arrivals of photons, 
which pertain to new fundamental physics, affecting the photon propagation. This may 
be due, for instance, to space-time foamy vacuum structures due to quantum gravitational 
fluctuations [9] that lead to modified dispersion relations for photons [101 EH E2] . If true, this 
would be a clear departure from the Lorentz invariant energy (E)-momentum (p) relations of 
Special Relativity, E = \p\c. 
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Figure 1. Comparison of data on delays At in the the arrival times of energetic gamma rays 
from various astrophysical sources with models in which the velocity of light is reduced by an 
amount linear in the photon energy. The graph plots on a logarithmic scale the quantity At/E 
and a function of the red-shift, K(z), which is essentially the distance of the source from the 
observation point. The data include two Active Galactic Nuclei (AGN), Mkn 501 [8j and PKS 
2155-304 [IB], and three Gamma Ray Bursts (GRB) observed by the Fermi satellite, 090510 [4J, 
09092B [5] and 080916c [3]. 



If we plot the observed time delays At versus the energy E of the photon, including 
experimental errors, we obtain [13] the diagram of fig. [TJ We observe from the figure that 
with the exception of GRB 090510 , observed by the FERMI satellite [1], the other four points 
can be fitted by a straight line: 

At/ El = (0.43 ± 0.19) * iTWs/GeV, Jf W S j* ^fAgf^r , (D 

assuming an expanding Universe within the framework of the standard Cosmological-constant- 
Cold-Dark-Matter (ACDM) model. The function K(z) expresses the effect of an expanding 
Universe onto the linear relation At oc E, and includes (i) a time dilation factor |12| (1 + z) and 
(ii) the redshifting [11] of the photon energy which implies that the observed energy of a photon 
with initial energy E is reduced to -E bs = Eq/(1 + z). Taking into account that in a Robertson- 
Walker expanding Universe, we assume for our analysis, the infinitesimal time interval dt is 
related to the Hubble rate H(z), where z is the red-shift, via: dt = — dz/[(l + z)H(z)], we obtain 
a total delay in the arrival times of photons with energy difference AE: 

(At) u - AE H" 1 P {i+z)dz ,<y, 

(Atj obs - Mqgi ±i Jo v / nA+nm( i +2) 3 W 



where Ho is the (current-era) Hubble expansion rate, which is fitted to the data, ([I]). As already 
mentioned, we have assumed for concreteness the ACDM standard model of cosmology, with 
= ^j^" representing the present-epoch energy densities (pj) of matter (including dark matter), 

3H 2 

f2 m , and dark vacuum energy, S7a, in units of the critical density p c = 87rG ° of the Universe 
(Gn is the Newton's gravitational constant). The current astrophysical measurements of the 
acceleration of the Universe are all consistent with a non zero Cosmological-Constant Universe 
with Cold-Dark-Matter (ACDM Model), with A ~ 73% and n m ~ 27%. The quantity M QG i 
appearing in (J2J) may be viewed as a phenomenological parameter at this stage, with units of 
energy. In concrete models of space-time foam, where propagation effects due to QG are assumed 
to be the dominant cause of the delay, this scale will be identified with the Quantum Gravity 
scale, where space-time foam effects are expected to set in. 

The best fit [13J, based on (TjQ), leads to the following result for the scale, Mqqi = 
(0.98^o 30) x 10 18 GeV. This result is remarkable, because it provides the first hint that such 
delays might be related with string theory effects, given that the order of the scale is that of the 
conventional string mass scale |14j . 

However, the GRB 090510 data, provided that the latter can be trusted 0, do not fit to this 
linear energy scheme. In fact the time delays and energies pertaining to the GRB 090510 [3] are 
such that the pertinent energy scale Mqqi in ([I]), that would reproduce this point in the graph 
of fig. CD is found to be at least 1.2M P , with M P ~ 1.2. x 10 19 GeV the Planck mass. 

The immediate reaction to this result, of course, is [1] to exclude the QG as being responsible 
for the induced time delays, since one expects, on accounts of naturalness, that any QG effect 
should not appear at scales higher than the Planck scale. And, indeed, this may well be the 
case, since, as we have already mentioned, astrophysical mechanisms at the source, which are 
in general different between GRBs and AGNs, may be responsible for the observed delayed 
arrivals of the more energetic photons. Nevertheless, as we shall discuss below, the pattern of the 
observed photon delays fits [T3] a string model of quantum-gravity-induced refractive index, with 
the pertinent quantum gravity energy scale being essentially the same as that inferred from the 
MAGIC observations (of order 10 18 GeV). Viewed lower bound, this scale is also compatible 
with that obtained from other Gamma- Ray data of the H.E.S.S. Collaboration [lb\ [Toj. In 
particular, we shall argue that a stringy model of space-time foam ("D-foam"), proposed in 
\17\ [T51 [T9] , which is based on scenarios involving brane worlds punctured by localized space- 
time brany defects (termed D- "particles" ) , can provide consistent fits to the current data on 
delayed arrival times of photons from MAGIC and FERMI Telescopes. 

However, in order to accommodate the findings of MAGIC with the recently observed time 
delays of photons from the extremely short Gamma-Ray burst GRB 090510 [4], one needs 
inhomogeneous densities of D-particle defects in the foam [20]. In this sense, astrophysical 
observations of such intense, short Bursts, are quite essential in falsifying models. The important 
feature of the string model, which allows for a consistent fit to all available data, including 
the GRB 090510, is that the foam there consists of real and not virtual space-time defects. 
Hence, their density may be non-uniform, varying with the redshift, i.e. be different at various 
epochs of the Universe. As we shall discuss below [20], the effective QG scale that enters ([2]) 
is inversely proportional to this density. A depletion or reduction, therefore, of the number of 
defects encountered by the photons at the redshift region of the GRB 090510, z = 0.9 (which 
by the way is the regime where one expects a cosmic deceleration to acceleration transition for 
the Universe) may well be responsible for a reconciliation of the observed delays for the GRB 
090510 with those of MAGIC, in the sense that they are both due to the (inhomogeneous) stringy 

1 Indeed, there are uncertainties in this measurement concerning, for instance, the precise emission time of 
photons due to the yet uncertain duration of pre-cursors to burst activities and other such issues, which need to 
be confirmed by other measurements of similarly short bursts, that are presently lacking. 



foam. From this point of view, the very few available data on observed time delays available 
today (fig. [1]) are not sufficient to exclude the model. One needs a significant improvement 
on statistics of relevant measurements at various redshifts and various directions (in order 
to exclude anisotropic foam situations), before definite conclusions on the falsification of the 
model are reached. Complementary tests from the Cosmology of these models may be helpful, 
although for inhomogeneous foam situations, studies in the early universe epoch cannot be used 
to constrain the present-era density of foam, as the two could be different. 



2. Lorentz- Violating Stringy Space-time Foam model: D-foam 

One of the cornerstones of Modern Physics is Einstein's theory of Special Relativity (SR), 
which is based on the assumption that the speed of light in vacuo c is an invariant under all 
observers. In fact, this implies invariance of the physical laws under the Lorentz transformations 
in flat space times, and the role of c as a universal limiting velocity for all particle species. The 
generalization (by Einstein) of SR to include curved space times, that is the theory of General 
Relativity (GR), encompasses SR locally in the sense of the strong form of the equivalence 
principle. According to it, at every space-time point, in an arbitrary gravitational field, it is 
possible to choose a locally inertia! ('free- float') coordinate frame, such that within a sufficiently 
small region of space and time around the point in question, the laws of Nature are described by 
special relativity, i.e. are of the same form as in unaccelerated Cartesian coordinate frames 
in the absence of Gravitation. In other words, locally one can always make a coordinate 
transformation such that the space time looks flat. This is not true globally, of course, and 
this is why GR is a more general theory to describe gravitation. The equivalence principle relies 
on another fundamental invariance of GR, that of general coordinate, that is the invariance of 
the gravitational action under arbitrary changes of coordinates. This allows GR to be expressed 
in a generally covariant form. 

In such a locally Lorentz-invariant vacuum, the photon dispersion relation, that is a local 
in space-time relation between the photon's four-wavevector components k^ = (uj,k) (where u 
denotes the frequency, and k the momentum) reads in a covariant notation: 

Wr}^ = (3) 

where repeated indices /x, v = 0,1,... 3, with referring to temporal components, denote 
summation and rj^ denotes the Minkowski space-time metric, with components r/oo = — 1, r/oi = 
Vio = 0) Vij = Vji = <% , * = 1,2,3 with Sij the Kronecker delta symbol. 

The above relation ([3|) implies the equality of all three kinds of photon velocities in vacuo 
that stem from its wave nature (due to the particle- wave duality relation): 



phase : v ph = g = ^ = c 

group : Ugr = || = jj-J-y = c , n gr {uj) = n(u) + uj^ 1 

front : v {Tont = c/n(oo) = c (4) 



since the phase and group refractive indices of the trivial vacuum equal unity n(uj) = n gr (w) = 1. 
For brevity we shall work from now on in units where % = c = 1. 

GR is a classical field theory. In a quantum theory of space-time, i.e in a Quantum-Gravity 
(QG) model, there is a priori no fundamental reason for the strong equivalence principle to 
hold, and in this sense, local Lorentz invariance might not be respected. This is the situation, 
in fact, that characterizes certain QG models, which predict a space-time "foamy" structure at 
microscopic (Planck) scales [9|. The latter may not be Lorentz invariant. The foamy structures 
could include a variety of topologically non-trivial configurations at such microscopic scales, 



ranging from virtual black holes to space-time non-commutativity and (real) stringy defects, as 
is the case of our stringy D-particle foam models [17} [TBI 119] . which we now proceed to review. 




Figure 2. A string theory model of D-particle "foam" . The model consists of appropriate stucks 
of parallel D(irichlet)branes, some of which are moving in a higher-dimensional Bulk space time, 
punctured by point-like D-brane defects (D-particles). The observable Universe is represented 
by one of such moving branes, compactified to three spatial dimensions (D3-brane). As the 
D3-brane world moves, D-particles from the bulk cross it and, thus, to an observer on the brane, 
they appear as "flashing on and off" space-time foam defects ("D-particle foam"). Photons, 
represented by open strings with their ends attached on the D3 brane, interact with these 
defects via capture/recoil, and this leads to non-trivial refractive indices. The effect is therefore 
"classical" from the bulk space time point of view, but appears as an effective "quantum foam" 
from the D3-brane observer effective viewpoint. 

The model of [T7] constitutes an attempt to construct a brane/string- inspired model of space 
time foam which could have realistic cosmological properties. For this purpose we exploited the 
modern approach to string theory [J3], involving membrane hypersurfaces (D(irichlet)-branes). 
Such structures are responsible for reconciliating (often via duality symmetries) certain string 
theories (like type I), which before were discarded as physically uninteresting, with Standard- 
Model phenomenology in the low-energy limit. 

In particular, we considered (c./. figure [2|) a ten-dimensional bulk bounded by two 
eight-dimensional orientifold planes, which contains two stacks of eight-dimensional branes, 
compactified to three spatial dimensions. Owing to special reflective properties, the orbifolds 
act as boundaries of the ninth-dimension. The bulk space is punctured by point-like DO-branes 
(D-particles), which are allowed in type IA string theory (a T-dual of type I strings |21j ) we 



consider in [T7] and here0. These are massive objects in string theory |14| . with masses M s /g s , 
where M s is the string mass scale (playing the role of the quantum gravity scale in string theory) , 
and g s < 1 is the string coupling, assumed weak for our purposes. These objects are viewed 
as space-time defects, analogous, e.g. to cosmic strings, but these are point-like and electrically 
neutral. I have to stress at this point that, according to modern ideas in string theory [II], the 
scale M s is in general different from the four-dimensional Planck-mass scale Mp = 1.2 x 10 19 
GeV/c 2 , and in fact it is a free parameter in string theory to be constrained by experiment. The 
energy scale M s c 2 can be as low as a few TeV; it cannot be lower than this, though, since if this 
were the case we should have already seen fundamental string structures experimentally. 

Supersymmetry dictates the number of D8-branes in each stack in the model, but does not 
restrict the number of DO-branes in the bulk. For definiteness, we restrict our attention for now 
to the type-IA model, in which the bulk space is restricted to a finite range by two appropriate 
stacks of D8-branes, each stack being supplemented by an appropriate orientifold eight-plane 
with specific reflecting properties, so that the bulk space-time is effectively compactified to a 
finite region, as illustrated in Fig. [2j We then postulate that two of the D8-branes have been 
detached from their respective stacks, and are propagating in the bulk. As discussed above, 
the bulk region is punctured by DO-branes (D-particles), whose density may be inhomogeneous. 
When there are no relative motions of the D3-branes or D-particles, it was shown in [IT] that 
the ground-state energy vanishes, as decreed by the supersymmetries of the configuration. Thus, 
such static configurations constitute appropriate ground states of string/brane theory. On 
the other hand, relative motions of the branes break target-space supersymmetry explicitly, 
contributing to the dark energy density. 

Our model assumes a collision between two branes from the original stack of branes (c./. 
figure [2]) at an early epoch of the Universe, resulting in an initial cosmically catastrophic 
Big-Bang type event in such non-equilibrium cosmologies [22]. After the collision, the branes 
bounce back. It is natural to assume that, during the current (late) era of the Universe, the 
D3-brane representing our Universe is moving slowly towards the stack of branes from which 
they emanated and the configuration is evolving adiabatically. Hence populations of bulk D- 
particles cross the brane worlds and interact with the stringy matter on them. To an observer 
on the brane the space-time defects will appear to be "flashing" on and off. The model we 
are using involves eight-dimensional branes and so requires an appropriate compactification 
scheme to three spatial dimensions e.g. by using manifolds with non-trivial fluxes (unrelated 
to real magnetic fields). Different coupling of fermions and bosons to such external fields 
breaks target space supersymmetry in a way independent of that induced by brane motion, 
which could be the dominant one in phenomenologically realistic models. The consequent 
induced mass splitting |23l |2"3] between fermionic and bosonic excitations on the brane world is 
proportional to the intensity of the flux field (a string generalization of the well-known Zeeman 
effect of ordinary quantum mechanics, whereby the presence of an external field leads to energy 
splittings, which are however different between (charged) fermions and bosons). In this way 
one may obtain phenomenologically realistic mass splittings in the excitation spectrum (at TeV 
or higher energy scales) owing to supersymmetry obstruction rather than spontaneous breaking 
(this terminology, which is due to E. Witten [25], means that, although the ground state could 
still be characterised by zero vacuum energy, the masses of fermion and boson excitations differ 
and thus supersymmetry is broken at the level of the excitation-spectrum). We also mention 
that the plausible assumption of a population concentration of massive D-particle defects in the 
haloes of galaxies can lead to modified galactic dynamics [26J . 

One can calculate the vacuum energy induced on the brane world in such an adiabatic 
situation by considering its interaction with the D-particles as well as the other branes in the 



2 One can extend the construction to phenomenologically realistic models of type IIB strings [19] . 



construction. This calculation was presented in [17\ I20j. where we refer the interested reader 
for further details. Here we mention only the results relevant for the present discussion. The 
important point to notice here is the fact that the forces exerted by the various structures, such 
as D-particles and other branes, on the brane world, representing the observable Universe in 
this kind of models, do not have a fixed sign. 

We concentrate first on D0-particle/D8-brane interactions in the type-IA model of |17j . 
During the late era of the Universe when the approximation of adiabatic motion is valid, we use 
a weak-string-coupling approximation g s <C 1. In such a case, the D-particle masses ~ M s /g s are 
large, i.e., these masses could be of the Planck size: M s /g s ~ Mp = 1.22- 10 19 GeV or higher. In 
the adiabatic approximation for the relative motion, these interactions may be represented by a 
string stretched between the DO-particle and the D8-brane, as shown in Fig. [21 The world-sheet 
amplitude of such a string yields the appropriate potential energy between the D-particle and 
the D-brane, which in turn determines the relevant contribution to the vacuum energy of the 
brane. As is well known p3], parallel relative motion does not generate any potential, and the 
only non-trivial contributions to the brane vacuum energy come from motion transverse to the 
D-brane. Neglecting a velocity-independent term in the D0-particle/D8-brane potential that is 
cancelled for a D8-brane in the presence of orientifold Og planes |14| we find |17| : 

V8S?D* = fOT r«Vc7, (5) 

VdTds = + i^7 for r»Vc7. (6) 



D0-D8 

where d < 1 is the relative velocity between the D-particle and the brane world, which is 
assumed to be non-relativistic. We note that the sign of the effective potential changes between 
short distances © and long distances ([6]). We also note that there is a minimum distance given 
by: 

r min ~ Vv~c7 , v < 1 , (7) 

which guarantees that ([5]) is less than r/a' , rendering the effective low-energy field theory well- 
defined. Below this minimum distance, the D0-particle/D8-brane string amplitude diverges 
when expanded in powers of (a') 2 v 2 /r 4 . When they are separated from a D-brane by a distance 
smaller than r m ; n , D-particles should be considered as lying on the D-brane world, and two 
D-branes separated by less than r m ; n should be considered as coincident. 

We now consider a configuration with a moving D8-brane located at distances Ri (t) from the 
orientifold end-planes, where Ri(t) + R2(t) = Ro the fixed extent of the ninth bulk dimension, 
and the 9-density of the D-particles in the bulk is denoted by n*(r): see Fig. The total 
D8- vacuum-energy density p 8 due to the relative motions is [17]: 



= - I n*{r) ^dr- n *(r) ^ dr + 

^ min 



* + h? W »*(»0£f7 dr + PO (8) 

where the origin of the r coordinate is placed on the 8-brane world and po combines the 
contributions to the vacuum energy density from inside the band — r m i n < r < r m i n , which 
include the brane tension. When the D8-brane is moving in a uniform bulk distribution of 
D-particles, we may set n* * (r) = no, a constant, and the dark energy density p^~i D0 on the 
D8-brane is also (approximately) constant for a long period of time: 

Ptoti D ° = -«o&t;(l -v) + n v 2 JE ^ 7 (R 1 (t) 2 + R 2 (t) 2 - 2a') + p . (9) 



3 This cancellation is crucial for obtaining an appropriate supersymmetric string ground state with zero ground- 
state energy. 



Because of the adiabatic motion of the D-brane, the time dependence of Ri{t) is weak, so that 
there is only a weak time dependence of the D-brane vacuum energy density: it is positive if 
po > 0, which can be arranged by considering branes with positive tension. 

However, one can also consider the possibility that the D-particle density n*(r) is 
inhomogeneous, perhaps because of some prior catastrophic cosmic collision, or some subsequent 
disturbance. If there is a region depleted by D-particles [20] - a D-void - the relative importance 
of the terms in §§§ may be changed. In such a case, the first term on the right-hand side of @ 
may become significantly smaller than the term proportional to Rf(t)/a'. As an illustration, 
consider for simplicity and concreteness a situation in which there are different densities of D- 
particles close to the D8-brane (ni oca i) and at long distances to the left and right (rMeft> n right)- 
In this case, one obtains from (JS]): 

Pto^ D ° - -rciocal^(l ~v)+ nMtV 2 ^ (R^t) 2 - a') + 

n right^l6lb {R2(tf ~ a') + p . (10) 

for the induced energy density on the D8-brane. The first term can be significantly smaller in 
magnitude than the corresponding term in the uniform case, if the local density of D-particles 
is suppressed. Overall, the D-particle-induced energy density on the D-brane world increases as 
the brane enters a region where the D-particle density is depleted. This could cause the onset of 
an accelerating phase in the expansion of the Universe. It is intriguing that the red-shift of GRB 
090510 [1] is in the ballpark of the redshift range where the expansion of the Universe apparently 
made a transition from deceleration to acceleration [27]. According to the above discussion, then, 
within our string foam model this may not be a coincidence [20 . The result (|10p was derived 
in an oversimplified case, where the possible effects of other branes and orientifolds were not 
taken into account. However, as we argued in [20J, the ideas emerging from this simple example 
persist in more realistic structures. 

The important point for our purposes in such models of space-time foam, is the fact that 
there is an induced refractive index in vacuo, as a consequence of the photon interactions with 
the D-particles that cross the moving brane Universe. It is the linear density of the D-particle 
defects n{z) encountered by a propagating photon that determines the amount of refraction. The 
density of D-particles crossing the D-brane world cannot be determined from first principles, 
and so may be regarded as a parameter in phenomenological models. The flux of D-particles 
is proportional also to the velocity v of the D8-brane in the bulk, if the relative motion of the 
population of D-particles is ignored. 

In order to make some phenomenological headway, we adopt some simplifying assumptions. 
For example, we may assume that between a redshift z < 1 and today (z = 0), the energy 
density has remained approximately constant, as suggested by the available cosmological data. 
As shown in [20J, this assumption implies certain restrictions for the bulk density w s hort( z ) of 
D-particle defects near the brane world 

"short (*) = "-short(O) - 277r io a ,5 jg f s jjZTj Jo H(z>) Z (l+z>) Where £s = Vc7- (11) 

We can make use of this result when we fit the available data from MAGIC and FERMI telescopes 
to this model. For instance, by requiring that this density falls from 0(1) at the redshifts 
z = 0.03, relevant to MAGIC observations, by at least two orders of magnitude as we approach 
the red-shift z = 0.9 of GRB 090510, we may constrain the relative motion of the D3-brane 
world. In fact, as discussed in [20], this can be achieved for relatively small velocities of the 
D3-brane world in the bulk, v < 10~ 4 c, consistent with the Early Universe Cosmology of the 
model [22]. 

Before closing this section we would like to make a comment regarding the modification of 
(thermal) Dark Matter (DM) relic abundances in the foam model. As discussed in [28j . the 



quantum fluctuations of the D-particles act as sources of particle production and thus affect the 
respective Boltzmann equation determining the relic abundance of DM particles. Additional 
modifications to this equation are due to the Finsler-like [29] character of the induced space- 
time metric during the interaction of neutral particles (like the DM ones) with the foam, due to 
its dependence on D-particle recoil velocities (and hence momentum transfer) (c.f. (|24p below). 
We assume foam types in which the Lorentz symmetry is conserved on the average, and is 
violated only through fluctuations, namely we assume statistical ensembles of D-particles such 
that <C Ui S>= 0, <C UiUj S>= j^jAfpiPj5ij, where p is some average momentum scale. One 
finds [28]: 
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where x = m/T, T is the temperature, the suffix / denotes the freeze-out point of the DM species, 
the suffix denotes present-day quantities, denote the Hubble-constant-free present-day 

thermal relic abundance of a weakly interacting heavy Dark Matter particle of mass m, in 
the presence of the foam, and the suffix "no source" indicates the corresponding quantity 
in the Standard (foam-free) Cosmology. Clearly the dominant correction terms are of order 

^(ELiA?)>0. 

The modifications are suppressed by the square of the string scale (actually the D-particle 
mass M s /g s ) and hence for relatively high string scales (much higher than TeV, of interest to us 
here) are small and do not lead to significant constraints on the density of foam. Indeed, in our 
case, from the MAGIC experiment, if we explain the observed delays as being due to the stringy 
foam alone, we obtain M s /g s ~ 10 19 GeV, while phenomenologically realistic DM candidates 
have masses m in the range of at most a few TeV (usually a few hundreds of GeV) |30j . 

It is important in this latter respect to make a comment regarding the nature of our foam: 
the D-particles in our approach \17\ [TBI [T9] are viewed as background configurations and not 
as excitations of the string vacuum. In other string/brane models, some authors have viewed 
the D-particles as localized excitations of the vacuum [31]. In those cases, the D-particles may 
be considered as dark matter candidates themselves and their density would be constrained 
by the cosmological observations on the DM sector to avoid overclosure of the Universe. As 
such they could not contribute to the refractive index. Indeed, for the latter property to occur, 
light must interact coherently with the D-particles, rather than scattering on them individually. 
Otherwise, there would not just be a time delay and thus an index of refraction, but the light 
would be incoherently deflected at arbitrarily large angles. Coherent scattering can only occur 
if the wavelength of light is much greater than the mean separation of the scatterers. In case 
the DO particles behave like dark matter [31], with masses near the Planck scale (to account for 
the MAGIC delays [8]), their number density must be less than 10 -20 m -3 to avoid overclosure 
of the Universe. Thus their mean separation would be greater than 100 m. The gamma rays 
that are of interest to us here have a much shorter wavelength (smaller than cm) and therefore 
cannot experience refraction due to these D-particles 0. In our D-foam model \\.7\ [T8] , where 
the DO-branes (or the compactified D3 branes around 3-cycles in the model of [19]) are viewed 
only as background defects, such an issue does not arise, and the density of the foam cannot 
be constrained by overclosure of the Universe issues. As we have discussed above, in this case 
the D-particles contribute to the dark energy sector (when in motion) consistently with current 
cosmological data. The cancellation between attractive (gravitational) and repulsive flux forces 
guarantees a supersymmetric vacuum, with zero vacuum energy, if no relative motion of foam 
occurs. 



4 At any rate, such super-heavy DM would have been washed out by inflation in any realistic cosmology, so the 
scenario of [3T] f° r D-particle excitations to play the role of DM pertains to much lighter D-particles in theories 
with low string mass scales. 



In the above picture, the presence of a ground-state D-particle defect would be reflected 
in the world-sheet boundary conditions for the open strings representing radiation or matter 
excitations. Thus, the only effects a recoiling D-particle would have can be summarized through 
its induced Finsler-like metric distortions and the associated time delays of the interacting 
neutral matter or radiation, which we now proceed to discuss. 



3. Time Delays in D-particle foam and Stringy Uncertainties 

In this section we discuss a possible origin of time delays induced in the arrival time of 
photons, emitted simultaneously from an astrophysical object, as a result of their propagation 
in the above-described D-particle space-time foam model. This comes by considering local 
interactions of photons, viewed as open string states in the model, with the D-particle defects. 
As we shall argue, the above model is in principle capable of reproducing photon arrival time 
delays proportional to the photon energies, of the kind observed in the MAGIC experiment [HE]. 
This microscopic phenomenon, which is essentially stringy and does not characterize local field 
theories, contributes to a sub-luminal non-trivial refractive index in vacuo, induced by the 
capture of photons or electrically neutral probes by the D-particle foam [18]. The capture 
process is described schematically in figure [3l An important feature of the model is that, on 
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Figure 3. Schematic view of the capture process of an open string state, representing a photon 
propagating on a D3-brane world, by a D-particle on this world. The intermediate string state, 
indicated by thick wave lines, which is created on capture of the end(s) of the photon by the 
D-particle, stretches between the D-particle and the brane world, oscillates in size between 
and a maximal length of order a'p°, where p° is the energy of the incident photon, and thus 
produces a series of outgoing photon waves, with attenuating amplitudes, constituting the re- 
emission process. The intermediate string state provides also the restoring force, necessary for 
keeping the D-particle roughly in its position after scattering. 

account of electric charge conservation, only electrically neutral excitations are subjected to 
capture by the D-particles. To the charged matter the D-particle foam looks transparent. This 
is because the capture process of fig. Oentails a splitting of the open string matter state. Charged 
excitations are characterized by an electric flux flowing across the string, and when the latter is 
cut in two pieces as a result of its capture by the D-particle defect, the flux should go somewhere 
because charge is conserved. The D-particle, being neutral, cannot support this conservation, 
and hence only electrically neutral excitations, such as photons, are subjected to this splitting 
and the associated delays. The reader should bear in mind of course that the D-particles carry 
other kinds of fluxes, unrelated to electromagnetism |14| . These are conserved separately, and it 
is for this reason that isolated D-particles cannot exist, but there must always be in the company 
of other D-branes, as in our model above, so that the relevant fluxes are carried by the stretched 
strings between the latter and the D-particles. 



For our purposes in this work it is also important to remark that the D-particles are treated 
as static when compared to photons. This is because the ends of the open string representing the 
photon move on the D3 brane world with the speed of light in (normal) vacuo, while the relative 
velocities of the D-particles with respect to the brane world (which propagates in the bulk 
space) are much lower than this. For instance, as discussed in [22], to reproduce cosmological 
observations in this model, in particular the spectrum of primordial density fluctuations, which 
are affected by the relative motion of the brane world, the speed of propagation of the D3-brane 
Universe should be smaller than 10 _4 c. 

When the end(s) of the open-string photon state are attached to the D-particle, there is 
an intermediate string state formed, stretched between the D-particle and the D3-brane, which 
absorbs the incident energy p° of the photon state to grow in size from zero to a maximum length 
L that is determined by the requirement of energy minimization as follows: one assumes [32] 
that the intermediate string state needs N oscillations to achieve its maximal length L, as in 
the standard string-string scattering case. Following the same logic, we first observe that energy 
conservation during the capture process of fig. [3J leads to relations of the form 

P° = £ + T (13) 

Minimizing the right-hand-side with respect to L and taking into account that the one end of 
the intermediate string state attached to the brane world moves with the speed of light in vacuo, 
c = 1 (in our units) , one arrives at time delays for the outgoing photon waves of order [18] : 

At ~ a'p° . (14) 

This delay is causal and can also be obtained by considering scattering of open string states off 
D-particle backgrounds □. One assumes that the incident open string state splits into two open 
string states upon interaction with the D-particle defect, which, after their scattering in the 
D-particle background, the presence of which is implemented through the appropriate Dirichlet 
world-sheet boundary conditions of the open strings, re-join to produce the outgoing state. 
The intermediate stage involves two open strings that scatter off each other. The (quantum) 
oscillations of the intermediate string state will produce a series of outgoing wave-packets, with 
attenuating amplitudes, which will correspond to the re-emission process of the photon after 
capture. The presence of the stretched string state, which carries the characteristic flux of the 
D-brane interactions, provides the restoring force, necessary to keep the D-particle in its position 
after scattering with the photon. 

The situation may be thought of as the stringy/brany analogue of the restoring force in 
situations in local field theories of photons propagating in media with non trivial refractive 
indices, as discussed by Feynman [33]. The role of the electrons in that case (represented 
as harmonic oscillators) is played here by the D-particle defects of the space-time. The 
stringy situation, however, is more complicated, since the D-particles have an infinite number 
of oscillatory excitations, represented by the various modes of open strings with their ends 
attached to them. Moreover, contrary to the conventional medium case, in the string model the 
refractive index is found proportional to the photon frequency, while the effective mass scale 
that suppresses the effect (fT4"|) is the quantum gravity (string) scale M s . The latter property can 
be understood qualitatively by the fact that the mass of these D-particle defects is of order |14] 
l/(g s Va*) = M s /g s , where g s is the string coupling (in units of h = c = 1). 

5 For type IA strings, which admit point-like D-particles, the situation is a straightforward extension of techniques 
applied to the open string-string scattering case of |32] , For type IIB strings, the situation is technically more 
involved, due to the compactified D3-branes, wrapped up around three cycles, that play the role of "D-particles" 
in the model [19] . but the results are qualitatively the same as in (|14[) . 



We remark at this stage that the above time delays are a direct consequence of the stringy 
uncertainty principles. Indeed, strings are characterized by two kinds of uncertainty relations: 
the phase-space Heisenberg uncertainty, modified by higher order terms in a' [34] as a result of 
the existence of the minimal string length £ s = yoi in target space-time, 

AXAP>h + a'(AP) 2 + ... (15) 

and the time-space uncertainty relation [35] 

AX At > a . (16) 

Since the momentum uncertainty AP < p , we have from (|15p . to leading order in a' (in units 

n = iy. 

^ x > > w ■ ^ 

In view of ()14|) . we then arrive at consistency with the space-time uncertainty (|16p . 

^ > ^ ~ £ (18) 

As in the conventional string theory photon-photon scattering [32], these delays are causal, i.e. 
consistent with the fact that signals never arrive before they occur. Hence they are additive for 
multiple scatterings of photons by the foam defects from emission till observation. As we shall 
discuss in section [U this provides the necessary amplification, so that the total delay of the more 
energetic photons can be |18] of the observed order in MAGIC and FERMI experiments. 

In addition to this leading refractive index effect (|14p . there are corrections induced by the 
recoil of the D-particle itself, which contribute to space-time distortions that we now proceed 
to discuss. From a world-sheet view point, the presence of D-particle recoil may be represented 
by adding to a fixed-point (conformal) a- model action, the following deformation [36, 37J: 

VrT = 2^7 Ef=i I dD dT Ui X°e (X*) d n X\ (19) 

where D in the sum denotes the appropriate number of spatial target-space dimensions. For a 
recoiling D-particle confined on a D3 brane, D = 3. 

There is a specific type of conformal algebra, termed logarithmic conformal algebra [38j , that 
the recoil operators satisfy |36} 137] , This algebra is the limiting case of world-sheet algebras 
that can still be classified by conformal blocks. The impulse operator Q(X°) is regularized so 
that the logarithmic conformal field theory algebra is respected [§. The conformal algebra is 
consistent with momentum conservation during recoil [36} [37] , which allows for the expression 
of the recoil velocity Uj in terms of momentum transfer during the scattering 

Ui = 9s p -^, (20) 

with being the D-particle "mass" and Ap = p\— P2 the associated momentum transfer of a 
string state during its scattering with the D-particle. 

We next note that one can write the boundary recoil/capture operator (|19p as a total 
derivative over the bulk of the world-sheet, by means of the two-dimensional version of Stokes 

6 This can be done by using the world-sheet scale, e~ 2 = ln(L/a) 2 , with a an Ultra- Violet scale and L the 
world-sheet area, as a regulator ggl [37]: O s (X") = -i ^-e luix0 . The quantity e -¥ 0+ at the end of the 
calculations. 



theorem. Omitting from now on the explicit summation over repeated i-index, which is 
understood to be over the spatial indices of the D3-brane world, we write then: 



£ (X )+X°5 £ (Xt 



d a X i (21) 



where 5 e (X°) is an e-regularized 5- function. For relatively large times after the impact at 
X° = (which we assume for our phenomenological purposes in this work), this is equivalent 
to a deformation describing an open string propagating in an antisymmetric -B^-background 
corresponding to an external constant in target-space "electric" field, 

B 0i ~ m , Bij = , (X° > 0) (22) 

where the X°5(X°) terms in the argument of the electric field yield vanishing contributions in 
the large time limit, and hence are ignored from now on. 

To discuss the space time effects of a recoiling D-particle on an open string state propagating 
on a D3 brane world, we should consider a a- model in the presence of the B-field (|22p . which 
leads to mixed-type boundary conditions for open strings on the boundary dT> of world-sheet 
surfaces with the topology of a disc. Absence of a recoil-velocity lij-field leads to the usual 
Neumann boundary conditions, while the limit where g^ u —> 0, with m ^ 0, leads to Dirichlet 
boundary conditions. 

In analogy with the standard string case in background electric fields [39, 32J, one obtains a 
non- commutative space-time if recoil of the D-particle is taken into account. This can be seen 
upon considering commutation relations among the coordinates of the first quantised a-model 
in the background (f22|) . As in the standard case of a constant electric background field, in 
the presence of a recoiling D-particle, the pertinent non commutativity is between time and 
the spatial coordinate along the direction of the recoil velocity field (for large times t after the 
impact at t = 0): 

\x\t]=ie 1 \ e 01 (=-e 10 ) = 9 = ± cT X, (23) 

where, for simplicity and concreteness, we assume recoil along the spatial X 1 direction. Thus, 
the induced non commutativity is consistent with the breaking of the Lorentz symmetry of the 
ground state by the D-particle recoil. The quantity U{ = ^ and u c = ^—7 is the Born-Infeld 
critical field. The space-time uncertainty relations (|23p are consistent with the corresponding 
space-time string uncertainty principle (|16p . 

Of crucial interest in our case is the form of the induced open-string effective target-space- 
time metric. The situation parallels that of open strings in external electric field backgrounds, 
discussed in refs. [32] 139] . Hence, the effective open-string metric, g°P en > electnc ; which is due 
to the presence of the recoil-velocity field u, whose direction breaks target-space Lorentz 
invariance, is obtained by extending appropriately the results of strings in constant electric 
field backgrounds [391 E2] to the background ([22]) : 

5 open,electric = ^ _ -2) ^ > ^,^ = 0,1 

^open.clectric = ^ = ^ Qther yalues _ ^ 

For concreteness and simplicity, we considered a frame of reference where the matter particle has 
momentum only across the spatial direction X , i.e. ^ p\ = p \\ u\ ,P2 = P3 = 0. Moreover, 
as in the standard case of strings in an electric field background, there is a modified effective 
string coupling 



(i-« 2 r < 25 > 



The fact that the metric in our recoil case depends on momentum transfer variables implies 
that D-particle recoil induces Finsler-type metrics [29], i.e. metric functions that depend on 
phase-space coordinates of the matter (photon) state. 

We now mention that the presence of the critical "background field" u c , at which both the 
metric (|24p and the effective coupling (j25[) vanish, is associated with the destabilization of the 
vacuum [40] when the field intensity approaches the critical value. Since in our D-particle foam 
case, the role of the 'electric' field is played by the recoil velocity of the D-particle defect, the 
critical field corresponds to the relativistic speed of light, in accordance with special relativistic 
kinematics, which is respected in string theory by construction. On account of (|20p . then, 
this implies an upper bound on the induced momentum transfer, and hence on the available 
momenta, for the effective field theory limit to be valid. Indeed, if we represent Ap in (|2(jp as 
a fraction of the incident momentum pi, Ap = rpi, r < 1, then the condition that the recoil 
velocity of the D-particle is below the speed of light in vacuo, as required by the underlying 
consistency of strings with the relativity principle, implies [U] 

Ap = r Pl < % => p x < % . (26) 

When the incident momentum approaches the order of this cutoff, the effective string coupling 
(|25p vanishes, while above that value the coupling becomes imaginary, indicating complete 
absorption of the string state by the D-particle. The space-time distortion due to recoil is so 
strong in such a case that there is no possibility of re-emergence of the string state, the defect 
behaves like a black hole, capturing permanently the string state. 

It is important to notice that in modern string theory the quantities M s ,g s are completely 
phenomenological. In fact, it is possible to construct phenomenologically realistic string/brane- 
Universe models, in the sense of being capable to incorporate the Standard Model particles at 
low energies, with string scales in such a way that M s /g s is significantly lower than the Planck 
scale. For instance, there are constructions [42J in the large extra dimension framework, for 
which M s /g s is of the order of 10 19 eV, i.e. of order of the conventional GZK cutoff in Lorentz 
invariant particle physics models. Thus, by embedding the D-particle foam to such models, 
one may have the appearance of a novel type Gretisen Zatsepin Kusmin (GZK) cutoff |43j . 
of similar order to the one obtained from conventional Lorentz invariance arguments, but of 
quite different origin: here it is the subluminal nature of the recoil velocity of the foam that 
sets the new upper bound in momentum transfer. This is not fine-tuning in our opinion, but 
indicates the appearance of a new upper bound in momenta, related implicity to the underlying 
Lorentz invariance of the string theory, which is broken spontaneously by the recoiling D-particle 
background. This point will become of importance later on, when we discuss constraints of our 
model coming from ultra- high-energy/infrared photon/photon scattering [44] . 

Before closing this section we make an important remark. The induced metric (|24p will affect 
the dispersion relations of the photon state: 

«^r) l/ n P eri ' electric — 

However, because the corrections on the recoil velocity U{ are quadratic, such modifications will 
be suppressed by the square of the D-particle mass scale. On the other hand, the presence of the 
D-particle recoil velocity will affect the induced time delays (|14|) by higher-order corrections of 
the form, as follows by direct analogy of our case with that of open strings in a constant electric 
field [32]: 

^^with D— foam recoil velocity — Oi i_^2 • 0^) 

i 

Thus, the D-particle recoil effects are quadratically suppressed by the D-particle mass scales, 
since (c.f. (|20p ) u\ oc g s Api/M s , with Api the relevant string-state momentum transfer. 



(27) 



However, the leading order delay effect (|28p . obtained formally by considering the limit of 
vanishing recoil velocity, is linearly suppressed by the string scale, and thus the induced time 
delays are disentangled in the string foam model from the modified Finsler dispersion relations 
()27p . This is important for the phenomenology of the model, as we shall discuss below. 
Lack of a Local Effective Field Theory Formalism in recoiling D-particle Foam: 
An important comment arises at this stage concerning the construction of a possible local 
effective field theory action in the case of D-particle foam, when recoil of the D-particle defect 
is considered. In view of the formal analogy of the problem with that of an external electric 
background field, one is tempted to apply the same considerations as those leading to the 
non-commutative (NC) effective actions in field theory, for instance for case of NC quantum 
electrodynamics of ref. [35]. Naively, one might think of writing down a Standard- Model- 
Extension type Lagrangian [36], including non-commutativity (and thus Lorentz- violating but 
CPT conserving) terms in the presence of D-foam background, with the a P parameter being 
replaced by 6° l = m = g s ^ i , the D-particle recoil velocity. 

However, this is not correct. As a result of the momentum-transfer dependence of the non- 
commutativity parameter in the D-foam case, it is not possible to write down the effective action 
in target space as a power series of local quantum operators. The momentum transfer is not 
represented by such operators when taking the Fourier transform. 

Thus writing down a local effective action for the recoiling D-particle in interaction with, 
say, a photon, is not possible. This has important consequences for phenomenology. Indeed, the 
time delays (|28|) . which are associated with the stringy uncertainty, are found proportional to the 
incident energy of the (split) photon state and are thus linearly suppressed by the string scale. 
The associated refractive index is therefore linearly suppressed, but this cannot be interpreted as 
an average propagation of photons in the context of some local effective action. As we have seen 
above, the associated anomalous photon dispersion terms in this case, induced by the Finsler 
metric (f!H|) . are quadratically suppressed by the string scale. 

Therefore, any analysis, such as those involving high energy cosmic rays [UJ using an effective 
linear dispersion relation, obtained from a local effective action, does not apply to our problem. 
Moreover, it is a generic feature of any local effective theory, that is a theory on flat space times 
involving higher derivatives operators, to yield birefringence, since the corresponding modified 
dispersion relations for photons reduce at the end of the day to a solution of the photon frequency 
as a function of the wave number k which is obtained from a quadratic equation. The latter 
admits two solutions with different propagation for the two photon polarizations. This is not 
the case of the uncertainty-related time delay (f28j) . which, as we have seen, is independent of 
the photon polarization [18] . 

4. Multiple Photon-D-particle Scatterings in the Foam and Total (observed) Time 
Delays 

The above-discussed time delays (fT3|) pertain to a single encounter of a photon with a D- 
particle. In case of a foam, with a linear density of defects n*/yo7, i.e. n* defects per string 
length, the overall delay encountered in the propagation of the photon from the source to 
observation, corresponding to a traversed distance D, is: 

Attotai = a/p°n*fy = ^n*D , (29) 

where p° denotes an average photon energy. When the Universe's expansion is taken into 
account, one has to consider the appropriate red-shift-z dependent stretching factors, which 
affect the measured delay in the propagation of two photons with different energies, as well as 
the Hubble expansion rate H{z). Specifically, in a Rob erston- Walker cosmology the delay due 
to any single scattering event is affected by: (i) a time dilation factor [12] (1 + z) and (ii) the 



redshifting [TT] of the photon energy which implies that the observed energy of a photon with 
initial energy E is reduced to E^bs = Eq/{1 + z). Thus, the observed delay associated with (fT4"j) 
is [IIIIE]: 

6t ohB = (1 + z)5t = (1 + z) 2 Vc7E ohs . (30) 

For a line density of D-particles n(z) at redshift z, we have n(z)d£ = n(z)dt defects per co- 
moving length, where dt denotes the infinitesimal Robertson- Walker time interval of a co-moving 
observer. Hence, the total delay of an energetic photon in a co-moving time interval dt is given 
by n(z)(l + z^C^/c/Eobsdt. The time interval dt is related to the Hubble rate H(z) in the 
standard way: dt = —dz/[(l + z)H(z)]. Thus, from (f3"0"j) we obtain a total delay in the arrival 
times of photons with energy difference AE, which has the form considered in 111] , namely 
it is proportional to AE and is suppressed linearly by the quantum gravity (string) scale, M s : 

(At) ob , = ^H -'ffn'( Z ) VnA '^,, (3D 

where z is the red-shift, Ho is the (current-era) Hubble expansion rate, and we have assumed for 
concreteness the ACDM standard model of cosmology. The reader should compare this relation 
with ([2]), based on naive considerations. The presence in (|3ip of the linear density of defects 
n*(z), which is in general red-shift dependent, is a crucial difference, and it is this feature that 
allows the string model to fit simultaneously the MAGIC result and the data from the short 
burst GRB 090510, as we shall discuss below. 

Notice in (|31 j) that the essentially stringy nature of the delay implies that the characteristic 
suppression scale is the string scale M s , which plays the role of the quantum gravity scale in 
this case. The scale M s is a free parameter in the modern version of string theory, and thus 
it can be constrained by experiment. As we have discussed in this article, the observations of 
delays of energetic (TeV) photons from AGN by the MAGIC telescope pQ can provide such an 
experimental way of constraining n* /M s in ([2]). For AE ~ 10 TeV, for instance, the delay ([2]) 
can lead to the observed one of order of minutes, provided M s /n* ~ 10 18 GeV (in natural units 
with c = 1) [8]. This implies natural values for both n* and M s , although it must be noted that 
n* is another free parameter of the bulk string cosmology model of [17J, considered here. In 
general, n*(z) is affected by the expansion of the Universe, as it is diluted by it, but also depends 
on the bulk model and the interactions among the D-particles themselves. For inhomogeneous 
foam situations, the dependence of n* on the redshift can be quite complicated [20J. For redshifts 
of relevance to the MAGIC experiments, z = 0.034 <C 1, one may ignore the z-dependence of n* 
to a good approximation. 

The total delay (|3ip may be thought of as implying [33J an effective subluminal refractive 
index n(E) of light propagating in this space time, since one may assume that the delay is 
equivalent to light being slowed down due to the medium effects. On account of the theoretical 
uncertainties in the source mechanism, however, the result of the AGN Mkn 501 observations of 
the MAGIC Telescope translate to upper bounds for the quantity n* /M s in ([SJ, which determines 
the strength of the anomalous photon dispersion in the string/D-particle foam model. 

In view of the above discussion, if the time delays observed by MAGIC can finally be 
attributed partly or wholly to this type of stringy space-time foam, and therefore to the stringy 
uncertainties, then the AGN Mkn 501, and other such celestial sources of very high energy 
photons, may be viewed as playing the role of Heisenberg microscopes and amplifiers for the 
stringy space-time foam effects. 

But how, then, can we incorporate the data from GRB 0905210 [4] in this model? It is 
clear that a fit with a linearly modified dispersion relation does not work in this case for the 
values of the Quantum Gravity scale M<^f IC ~ 10 18 GeV that fit the MAGIC, H.E.S.S. and 
the other FERMI data. Indeed, the observed short delays of this burst can be explained on the 



basis of linearly modified dispersion relations only if the quantum gravity scale is larger than 
Mqgi ^ 1.2M P . 

As we have already mentioned, this argument has been used in [4J in order to exclude all 
models of QG involving linearly suppressed time delays for photons, on account of naturalness. 
We would not agree with this statement. Leaving aside the fact that from a single measurement, 
with uncertainties on the precursor of the GRB, one cannot draw safe conclusions, we mention 
that what one calls a natural scale of QG is highly model dependent. For instance, as we have 
seen above, in the string foam model, the relevant "QG scale", that dictates the order of the 
foam-induced time delays of photons, is a complex function of the model's parameters and is 
not simply given by the string scale M s [20J. Indeed, as becomes clear from (|'29p (or ([2])) the 
relevant scale is not simply the D-particle mass, M s /g s , but a combination 

M QG -D-foam ~ g s n*(z) ( 32 ) 

involving the linear number density n*(z) of the foam defect, encountered by the photon during 
its propagation from the source till observation. This quantity depends on the bulk density of 
the D-particles, which in the model of |17| is a free parameter. 

Inhomogeneous bulks are perfectly consistent background configurations for our brane world 
scenario |20j. In order to match the photon delay data of Fig. Q] with the D-foam model, we 
need a reduction of the linear density of defects encountered by the photon by about two orders 
of magnitude in the region 0.2 < z < 1, whereas for z < 0.2 there must be, on average, one 
D-particle defect per unit string length l s . We therefore assume that our D-brane encountered a 
D-void when 0.2 < z < 1, in which there was a significant reduction in the bulk nine-dimensional 
density of defects. Such assumptions can be consistent with cosmological considerations on the 
dark sector of the model, as explained in |20j . In particular, an important parameter for the 
cosmology of the model is the propagation velocity v of the brane world in the bulk (c./. fig. [2]). 
As we have discussed previously (fTTj) this is related [20] to the density of the defects near the 
brane world, n short (z), and hence to the linear density of defects encountered by the photon on 
the brane world. 

To have a reduced density by two order of magnitude at redshifts z = 0.9, while having a 
density of defects of O(l) per string length at redshifts z < 0.1, one must consider the magnitude 
of v as well as the string scale i s . For instance, it follows from (|lip that for string energy scales 
of the order of TeV, i.e. string time scales £ s /c = 10~ 27 s, one must consider a brane velocity 
v < vl0 x 10~ n c, which is not implausible for a slowly moving D-brane at a late era of the 
Universe 0. This is compatible with the constraint on v obtained from inflation in [22], namely 
v 2 < 1.48 x lO^ 5 ^" 1 , where g s < 1 for the weak string coupling we assume here. On the other 
hand, if we assume a 9-volume Vg = (K£ s ) 9 : K ~ 10 3 and £ s ~ 10~ 17 /GeV, then the brane 
velocity v < 10~ 4 c. In our model, due to the friction induced on the D-brane by the bulk 
D-particles, one would expect that the late-epoch brane velocity should be much smaller than 
that during the inflationary era immediately following a D-brane collision [22] , 

We mention for completion, that the above considerations pertain strictly speaking to point- 
like D-particles, which are allowed only in type IA string constructions. However, as discussed 
in |20| . similar conclusions can be applied to type-IIB constructions of D-foam |19j . 

5. Other (astrophysical) constraints on quantum- gravity foam 

The sensitivity of the MAGIC (and FERMI) observations to Planck scale physics, at least for 
linearly suppressed modified dispersion relations, calls for an immediate comparison with other 
sensitive probes of non-trivial optical properties of QG medium. 

7 Much smaller velocities are required for small string scales that are comparable to the four-dimensional Planck 
length. 



Indeed, from the analysis of [8], there was no microscopic model dependence of the induced 
modifications of the photon dispersion relations, other than the sub-luminal nature of the 
induced refractive index and the associated absence of birefringence, that is the independence 
of the refractive index on the photon polarization. The latter feature avoids the otherwise very 
stringent constraints on the photon dispersion relation imposed by astrophysical observations, 
as we now come to discuss. 

We shall be very brief in our description of the complementary astrophysical tests on Lorentz 
invariance and quantum-gravity modified dispersion relations, to avoid large diversion from our 
main point of this review article which is string theory. 

There are three major classes of complementary astrophysical constraints, to be considered 
in any attempt to interpret the MAGIC, FERMI or more general 7-ray Astrophysics results in 
terms of quantum-gravity induced anomalies in photon dispersion. 

• Birefringence and strong constraints on QG -induced photon dispersion 

In certain models of quantum gravity, with modified dispersion relations, for instance the 
so-called loop-quantum gravity [17], the ground state breaks reflexion symmetry (parity) 
and this is one of the pre-requisites for a dependence of the induced refractive index on the 
photon polarization, i.e. birefringence. We remind the reader that in birefringent materials 
this is caused precisely by the existence of some kind of anisotropies in the material. The 
velocities of the two photon polarizations (denoted by ±) in such QG models may be 
parametrised by: 



where M P = 1.22 x 10 19 GeV is the Planck energy scale, and £ is a parameter following from 
the underlying theoretical model, which is related with the modifications of the pertinent 
dispersion relations for photons. The order of suppression of these effects is described by n 
which in the models of [47] assumed the value n = 1, but in general one could have higher 
order suppression. 

Vacuum QG birefringence should have showed up in optical measurements from remote 
astrophysical sources, in particular GRBs. Ultraviolet (UV) radiation measurements from 
distant galaxies [38] and UV/optical polarization measurements of light from Gamma Ray 
Bursters [49J rule out birefringence unless it is induced at a scale (way) beyond the Planck 
mass (for linear models, the lower bound on the QG scale in such models can exceed the 
Planck scale (~ 10 19 GeV) by as much as seven orders of magnitude). Indeed, in terms 
of the parameter n introduced above (c.f ([33])), for the case n = 1 of [47] . one finds from 
optical polarization observations that the absence of detectable birefringence effects imply 
the upper bound |£| < 2 x 1CP 7 , which is incompatible with the MAGIC observed delays, 
saturating from below the bounds imposed by the MAGIC experiment. 
At this point, we wish to mention that, using recent polarimetric observations of the Crab 
Nebula in the hard X-ray band by INTEGRAL [50] . the authors of [51] have demonstrated 
that the absence of vacuum birefringence effects constrains linearly suppressed Lorentz 
violation in quantum electrodynamics to the level |£| < 6 x 1CP 10 at 95% C.L., thereby 
tightening by about three orders of magnitude the above-mentioned constraint. 

• Synchrotron Radiation and further stringent constraints for electronic QG- 
induced anomalous dispersion in vacuo 

Another important experimental constraint on models with QG-induced anomalous 
dispersion relations comes from observations of synchrotron radiation from distant 
galaxies [52l ESI EU > such as Crab Nebula (c./. fig. S])- As well known, the magnetic fields 
at the core regions of galaxies curve the paths of (and thus accelerate) charged particles, 
in particularly electrons (which are stable and therefore appropriate for astrophysical 




(33) 



observations), and thus, on account of energy conservation this results in synchrotron 
radiation. 




Figure 4. The Crab Nebula is a supernova remnant, with a rotating neutron star (the Crab 
Pulsar) at its centre (middle image). Observations of synchrotron radiation from such celestial 
objects place very stringent constraints on quantum gravity models with anomalous dispersion 
relations for electrons (from the first of Ref. |53j). 

In standard electrodynamics [55], electrons in an external magnetic field H, follow helical 
orbits transverse to the direction of H. The so-accelerated electrons in a magnetic field 
emit synchrotron radiation with a spectrum that cuts off sharply at a frequency uj c (c.f. 
fig- HI right panel): 

LI _ 3 eH 1 foA\ 

where e is the electron charge, mo its mass, and = v± is the component of the velocity 
of the electron perpendicular to the direction of the magnetic field. The superfix LI in (|34p 
stresses that this formula is based on a LI approach, in which one calculates the electron 
trajectory in a given magnetic field H and the radiation produced by a given current, using 
the relativistic relation between energy and velocity. 

All these assumptions are affected by violations of Lorentz symmetry, such as those 
encountered in quantum-gravity space-time foam models, leading to modified dispersion 
relations of the form: 

u> 2 (k) = k 2 + ^^-, (35) 
E\p) = m2+p 2 + e e ^J, (36) 

for photons (|35|) and electrons (|36|) . where oj and k are the photon frequency and wave 
number, and E and p are the electron energy and momentum, with mo the electron (rest) 
mass. In the spirit of the MAGIC observation analysis above, we assume here linear (a = 1) 
or quadratic QG (a = 2) effects, characterized by parameters £ 7 and £ e , extracting the 
Planck mass scale Mp = 1.22 x 10 19 GeV. In fact one can do the analysis [SS] for a 
general a (single power) and attempt to extract limits on this parameter by matching with 
observations. 

A detailed analysis (52J [53] , including the modifications in the electron's trajectories due to 
space-time foam [53], yields: 

u? G oc ^ 1 



( i +> /2=i7^ ) i/ a (^ + (a + l)(^ 



M = Mp/\C e \ , 7? = 1 - (£/.M) 



(37) 



where u^ 1 is given in (|34p and the superscript "QG" indicates that the QG-modified 
dispersion relations (|36p are used. This function is plotted schematically (for a = 1) in 
fig. H] (right panel). 

In |52l 153]. the above QG-modified dispersion relations have been tested using observations 
from Crab Nebula. It should be emphasized that the estimate of the end-point energy 
of the Crab synchrotron spectrum and of the magnetic field used above are indirect values 
based on the predictions of the Synchrotron Self-Compton (SSC) model of very-high-energy 
emission from Crab Nebula [SU] . In [S3] the choice of parameters used was the one that gives 
good agreement between the experimental data on high-energy emission and the predictions 
of the SSC model [56[ 157] . Estimating the magnetic field of Crab Nebula in the region 
160 x 10~ 6 Gauss < H < 260 x 10~ 6 Gauss, and requiring |£ e | < 1 (which thus sets the 
quantum gravity scale as at least Mp) one obtains the following bounds on the exponent a 
of the dispersion relations ([36]) [53] : 

a>a c : 1.72 < a c < 1.74 (38) 

These results imply already a sensitivity to quadratic QG corrections with Planck mass 
suppression Mp. 

However, for photons there are no strong constraints on £ 7 coming from synchrotron 
radiation studies, unless in cases where QG models entail birefringence [53], where, as 
we discussed above, strong constraints on the photon dispersion are expected at any rate 
from optical measurements on GRBs. In this sense, the result (|38p excludes the possibility 
that the MAGIC observations leading to a four-minute delay of the most energetic photons 
are due to a quantum foam that acts universally among photons and electrons. However, 
the synchrotron radiation measurements cannot exclude anomalous photon dispersion with 
linear Planck-mass suppression, leading to a saturation of the lower bound imposed by the 
MAGIC experiment [8], in models where the foam is transparent to electrons, as in the 
string foam case. 

• Strong constraints from Ultra- high- energy Cosmic photon annihilation 

Further strong constraints on generic modified dispersion relations for photons, like the 
ones used in the aforementioned QG-interpretation of the MAGIC results [8], comes from 
processes of scattering of ultra-high-energy photons, with energies above 10 19 eV, off very- 
low energy cosmic photons, such as the ones of the cosmic microwave background (CMB) 
radiation that populates the Universe today, as a remnant from the Big-Bang epoch. In 
[44] it has been argued that the non-observation of such ultra-high energy (UHE) photons 
places very strong constraints on the parameters governing the modification of the photon 
dispersion relations, that are several orders of magnitude smaller than the values required 
to reproduce the MAGIC time delays, should the effect be attributed predominantly to 
photon propagation in a QG dispersive medium. 

The main argument relies on the fact that an ultra-high-energy photon would interact with 
a low-energy ("infrared") photon of the CMB background (with energies in the eV range) 
to produce electron prositron pairs, according to the reaction: 

7UHE + 7CMB e + e~ . (39) 

The basic assumption in the analysis is the strict energy and momentum conservation in 
the above reaction, despite the modified dispersion relations for the photons. Such an 
assumption stems from the validity of a local-effective-lagrangian description of QG foam 
effects on particles with energies much lower than the QG energy scale (assumed close to 



Planck scale Mp\ = 10 19 GeV). In this formalism, one can represent effectively the foam 
dispersive effects by higher-derivative local operators in a flat- space-time Lagrangian. The 
upshot of this is the modification of the pertinent equations of motion for the photon field 
(which in a Lorentz-invariant theory would be the ordinary Maxwell equations) by higher- 
derivative terms, suppressed by some power of the QG mass scale. 

One considers the modified dispersion relations (|35p . (|36p. which in the notation of [44] . 
taking explicit account of the various polarizations and helicities, can be written as: 

ul = k 2 + £k 2 (^) n , u 2 = k 2 , 

El ± =p 2 e + m 2 e + rt ± p 2 e (^) n (40) 

with (co,k) the four-momenta for photons, and (E,p e ) the corresponding four-momentum 
vectors for electrons; the suffix b indicates a low-energy CMB photon, whose dispersion 
relations are assumed approximately the normal ones, as any QG correction is negligible due 
to the low values of energy and momenta. The +(-) signs indicate left(right) polarizations 
(photons) or helicities (electrons). Positive (negative) £ indicate subluminal (superluminal) 
refractive indices. Upon the assumption of energy-momentum conservation in the process, 
one arrives at kinematic equations for the threshold of the reaction (|39p , that is the minimum 
energy of the high-energy photon required to produce the electron-positron pairs. 
For the linear- or quadratic- suppression case, for which n = 1,2 respectively (|40p . one 
finds that, for the relevant subluminal photon refractive indices corresponding to the 
saturation of the lower bound on the QG scale Mqgi ~ 10 18 GeV, the threshold for pair 
production disappears for ultra- high-energy photons, and hence such photons should have 
been observed. The non-observation of such photons implies constraints for the relevant 
parameters £, rj which are stronger by several orders of magnitude than the bounds on the 
QG scale inferred from the MAGIC observations. 

From the analysis of [H] one concludes that in the case of linear Planck-mass suppression 
of the sub-luminal QG-induced modified dispersion relations for photons, of interest for 
the QG-foam interpretation of the MAGIC results [SJ, parameters with size £i > 10 
are ruled out. This exceeds the sensitivity of the MAGIC experiment to such Lorentz- 
symmetry violating effects by fifteen orders of magnitude ! Similar strong constraints are 
also obtained from the non observations of photon decay (7 —> e + e~), a process which, if 
there are modified dispersion relations, is in general allowed [44] . 



6. The String Foam Models Evade the above Constraints 

From the above discussion it becomes clear that any model of refraction in space-time foam 
that exhibits effects at the level of the MAGIC experiment sensitivity [8] should be characterised 
by the following specific properties: 

• (i) photons are stable (i.e. do not decay) but should exhibit a modified subluminal refractive 
index with Lorentz- violating corrections that grow linearly with E/(Mqq^c 2 ), where Mqq-j 
is close to the Planck scale, 

• (ii) the medium should not refract electrons, so as to avoid the synchrotron-radiation 
constraints [53] , and 

• (iii) the coupling of the photons to the medium must be independent of photon polarization, 
so as not to have birefringence, thus avoiding the pertinent stringent constraints [48} l4"9l 

• (iv) The formalism of local effective lagrangians should break down, in the sense that there 
are quantum fluctuations in the total energy in particle interactions, due to the presence 



of a quantum gravitational 'environment', such that stringent constraints, which otherwise 
would have been imposed from the non-observation of ultra-high energy photons (%uj > 10 19 
eV), are evaded. 

The string-foam models (both type IA and type IIB) are characterized by all these 
properties [T71 CEB] , and thus avoid the stringent constraints. The absence of birefringence 
and the transparency of foam to electrons or charged probes, for reasons of charge conservation, 
make the models surviving the stringent constraints from synchrotron radiation from distant 
Nebulae [53 GUI- 

Moreover, the analysis in [33] is based on exact energy momentum conservation in the process 
(|39p . stemming from the assumption of the local-effective lagrangian formalism for QG-foam. As 
we discussed in |58}l41j. and mentioned in section [31 however, such a formalism is not applicable 
in the case of the recoiling D-particle space-time foam model, where the fluctuations of space- 
time or other defects of gravitational nature paly the role of an external environment, resulting 
in energy fluctuations in the reaction (|39|) . The presence of such fluctuations does affect the 
relevant energy-threshold equations, for the reaction to occur, which stem from kinematics, in 
such a way that the above stringent limits are no longer valid. In particular, we have seen that 
the pertinent anomalous dispersion terms are induced by the Finsler-like metric (|24p . which in 
turn arises by the distortion of space time due to the recoil of the D-particle space-time defect 
during its interaction with the photon. This metric is quadratically suppressed by the string 
scale, and hence the so-induced modified dispersion relations (|27p contain anomalous terms 
quadratically suppressed by the QG scale, a case for which the constraints coming from ultra 
high energy cosmic rays are currently weak, but of course the situation can change in future 
experiments [44] . 

It should be stressed again that in the string model the modified dispersion relations are 
disentangled from the time delays (|28p of the more energetic photons, which are linearly 
suppressed by the string scale. These delays are associated with the string uncertainty 
principle [35J and thus are not represented within the local effective field theory framework. We 
have also seen that this disentanglement allows for a consistent interpretation of the MAGIC 
delays with those of the short GRB 090510 observed by FERMI within the string D-foam 
framework, provided of course the foam is inhomogeneous. In addition, as discussed above, the 
form of the effective string coupling ([25]) , and the induced metric ([24")) , imply an upper bound 
in the photon momentum transfer (|26p . which characterizes the D-foam model as a result of 
the sub-luminal nature of the D-particle recoil velocity. These considerations imply that very 
high energy cosmic rays, when interacting with the D-foam, will have an extremely suppressed 
interaction rate (the string amplitudes are vanishing when the recoil velocity approaches the 
speed of light), thus providing additional reasons [41j for evading the strong constraints of |44j . 
There are also phenomenologically realistic brane models, with large extra dimensions [32]) for 
which the ratio M s /g s [i.e. in our context the D-particle mass) may be of the order of the 
conventional GZK cutoff [43J, 10 20 eV. In these models, in view of (|26p . there would be no 
photons with energies higher than this that could not be completely absorbed when interacting 
with the D-particle foam. This, therefore, provides an explanation for the absence of such high- 
energy photons, in accordance with observations, and enables the model to evade the stringent 
constraints coming from ultra-high-energy cosmic rays. 

7. Conclusions and Outlook 

In this review we have discussed a stringy version of a (Lorentz-Invariance- Violating) space- 
time foam model, which appears as a candidate theory for an explanation of the delays of the 
more energetic photons from celestial sources, as observed by MAGIC and FERMI Telescopes, 
in agreement with all the other current astrophysical tests of Lorentz Violation. Although such 



delays could be due to conventional astrophysics at the source, nevertheless the latter is not well 
understood at present, and in fact there is no consensus among the relevant communities. 

The fact that string theory (or better its modern version involving D-brane defects) is 
capable of explaining potentially the observed photon delays, in agreement with all the other 
astrophysical data currently available, is at least amusing. This means, if nothing else, that 
at least such tests provide a framework for experimentally testing/falsifying some models of 
string theory (entailing Lorentz Violation) at present or in the foreseeable future. The key 
point in the approach is the existence of space-time defects in the ground state of the model, 
whose topologically non-trivial interactions with the string states, via string-stretching during 
the capture process (c.f. fig. E]), are mainly responsible for the observed delays. The latter are 
found proportional to the incident photon energy. The situation may be graphically represented 
by the case of two boats, of different lengths, in rough seas: the waves in the sea represent 
the foam, the "medium" over which the two boats propagate; the longer boat, representing 
the longer wave-length (lower energy) photons, will cross the waves faster than the smaller 
boat, corresponding to the higher energy (shorter wavelength) photons, which will be relatively 
delayed by having to climb up and down the waves. 

The peculiarity of the D-particle foam in being transparent to charged particles (as a result of 
electric charge conservation requirements), evades the stringent constraints on linear Planck scale 
suppression refractive indices that would otherwise have been induced by electron synchrotron 
radiation studies from Crab Nebula [52\ 154] . Moreover, the absence of birefringence avoids the 
similarly stringent constraints on such models that would have been imposed by galactic [S3] or 
extra-galactic measurements [48 \ H9j . Finally, it worths mentioning that the linear in energy time 
delays (fl~4"|) .([2j). when applied to neutrinos, can be flavour (i.e. neutrino-species) independent, 
thus avoiding [59] the stringent constraints that would be obtained from models of quantum 
gravity with flavour-dependent modifications of neutrino propagation and thus modifications in 
their oscillations [60J. 

We must stress once more, however, that despite the apparent "success" of the D-foam 
model in dealing with stringent astrophysical constraints so far, one cannot as yet draw any 
safe conclusions on the validity of the model in Nature. Many more high-energy astrophysical 
photon measurements are needed in order to disentangle source from possible propagation effects 
due to fundamental physics. If a statistically significant population of data on photons from 
cosmic sources is collected, exhibiting refractive indices varying linearly with the distance of the 
source [JJT], as well as the photon energy, then this would be a very strong confirmation of the 
D-particle foam model, for reasons explained above. However, it must be noted that GRB's, 
which are expected to lead to statistically significant data in the next few years, will produce 
photons much lower in energies than the flares observed in AGN, and this could be a drawback. 
At any rate there are attempts to claim that observations from FERMI will have sensitivity 
close to the Planck scale [5T] for such linear-suppression models soon. The case of GRB090510 
is a perfect example of how a single measurement of a short, intense high energy burst, can 
place stringent limits on Lorentz Violation at the Planck scale. However, caution should be 
exercised here as to what one means by "sensitivity at the Planck scale". As demonstrated 
above (c.f. (|32p). the effective quantum gravity scale is actually a complicated function of many 
microscopic parameters in the model. Nevertheless, if many observations at various redshift 
regimes on delayed arrivals of cosmic photons become available in the future, then we shall be 
able to make some definite conclusions regarding the order of magnitude of possible quantum 
gravity effects and thus falsify models, such as the D-particle space-time foam. 

We believe we are entering exciting times for tests of fundamental physics concepts and 
symmetries by means of high-energy astrophysics. Therefore it seems that this branch of physics, 
together with collider and particle physics, will provide a platform for exploring our world at 
increasingly smaller microscopic length scales, thereby offering experimental guidance in our 



quest for the elusive theory of Quantum Gravity. 
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